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Abstract 

The periodic standing wave method studies circular orbits of compact objects coupled to helically 
' symmetric standing wave gravitational fields. From this solution an approximation is extracted 

, for the strong field, slowly inspiralling motion of black holes and binary stars. Previous work on 

this model has dealt with nonlinear scalar models, and with linearized general relativity. Here we 
present the results of the method for the post-Minkowski (PM) approximation to general relativity, 
^ , the first step beyond linearized gravity. We compute the PM approximation in two ways: first, 

via the standard approach of computing linearized gravitational fields and constructing from them 
quadratic driving sources for second-order fields, and second, by solving the second-order equations 
as an "exact" nonlinear system. The results of these computations have two distinct applications: 
(i) The computational infrastructure for the "exact" PM solution will be directly applicable to full 
^ ■ general relativity, (ii) The results will allow us to begin supplying initial data to collaborators 

running general relativistic evolution codes. 
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I. BACKGROUND AND INTRODUCTION 



The inspiral of binary black holes is of great interest as a source of detectible gravitational waves and the computation 
of the waves from the inspiral has been the focus of much effort. Recent breakthroughs in numerical relativity P, Q 
hold the promise of computing the evolution of the last few orbits of inspiral of a binary pair of black holes. What 
remains is to find results for the epoch of inspiral earlier than the last few orbits, and to provide optimal initial data 
' for the evolution equations. 

QQ The Periodic Standing Wave (PSW) project is intended to fill this gap in a more-or-less efficient way. This method 

seeks a numerical solution for a pair of sources (black holes, neutron stars) in nondecaying circular orbits with 
, gravitational fields that are rigidly rotating, that is, fields that are helically symmetric. Because the universality of 
' gravitation will not permit outgoing waves and nondecaying orbits, the solution to be computed is that for standing 
^ ' waves. An approximation for slowly decaying orbits with outgoing radiation is then extracted from that numerical 
•""j ■ solution. 

This work has progressed through several stages. In the first staged, 13, S Hi 13 1 * nonlinear scalar fields model 
^ . was investigated, and numerical methods were developed to deal with the special mathematical features that would 
" " ' be common to all standing-wave, helically symmetric computations. These features include: (i) a mixed boundary 
value problem (regions of the domain in which the equations are hyperbolic and other regions in which they are 
elliptic); (ii) an iterative construction of nonlinear standing wave solutions; (iii) the extraction from the standing wave 
solution of an approximate outgoing wave solution; (iv) the effectiveness of Newton-Raphson methods to deal with 
the nonlinearities. In Ref. ^ standard finite-difference methods were used to explore the nonlinear scalar problem, 
but it was apparent that sufhcient resolution to achieve good convergence of the nonlinear iterations would involve 
a computationally intensive project, something we wanted to avoid. Reference Q introduced a new technique for 
greatly reducing the computational burden. That reference introduced "adapted coordinates" that were well suited 
to the geometry of the problem. Near each of the sources these coordinates approached spherical coordinates centered 
on the source; far from the sources, the coordinates approached standard spherical coordinates centered on the center 
of mass. In Ref. it was shown that with these coordinates good results could be computed by keeping only a very 
small number of multipoles, typically just the monopole and quadrupole moments. 

With the mathematical and computational methods for scalar fiels under some control we turned to linearized 
gravity in the harmonic gauge Q. The goal in that work was to describe linearized gravity with convenient "helical 
scalars" (functions only of coordinates corotating with the helical Killing congruence). We presented a formalism that 
was remarkably simple. Metric perturbations in the harmonic gauge were described as fields in a Minkowski spacetime, 
three complex fields vj/^^i)^ $(22)^ ^^^^ ^^^^ ^^^^ gg^jg ^(nn)^ §(nO)^ ^^(oo)^ §(20)^ ^j^j^ ^j^jg description, each of 
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the equations of hnearized general relativity, for each of fields, was found to have the form 

l(vI/('^'^)) = (1) 

outside sources. For the four real fields v^t'^"', the operator L is simply the second-order operator of linear scalar 
theory. For the three complex fields, L has extra terms. One of the extra terms turns out to be imaginary, so that 
the real and imaginary parts of the complex fields v[;("i)^ xi/(2i)^ v[;(22) g^j.g coupled by the equations. But this is the 
only coupling that occurs in this infrastructure for helically symmetric linearized general relativity. The presence of 
the sources is introduced in Ref. 8] through inner boundary conditions on a small, approximately spherically surface 
in the source region. 

Here we take the gravitational problem beyond linear theory, by truncating at second order a nonlinear expansion of 
the vacuum general relativistic field equations. We take two distinctly different views of the post-Minkowski problem 
that results. The first, to be called the "perturbative post-Minkowski" (PPM) problem, is the standard approach to 
a sequence of perturbation orders. In this approach, the field equations are written in the form 

= (2) 
^(f''^'''') = S'''^") ^^f^"', ^^''''^y (3) 

The notation here indicates that one is first to solve the linearized equations for the first-order perturbative fields 
One then constructs effective sources for the second-order fields . These effective sources are quadratic 

in the linear fields . The inner boundary conditions (representing the physical sources) must only be correct to 

linear order in solving Eq. ([2|) for ^'(^'') . In Eq. ^ only the second-order part of the boundary value must be used 
for . 
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The other approach to the second-order post-Minkowskian solution is to replace Eqs. ^ and ([3]), by a single 
equation 
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In this formulation, the field equations of general relativity are truncated at second order in the field strengths, and 
the resulting system, quadratic in the fields, is treated as a nonlinear field theory and solved as such. In this "exact 
post-Minkowski" (exact PM) approach, in contrast to that of Eq. ([3]), there is no a priori division of v[;('»f') into first- 
and second-order parts, and the boundary conditions for vj/l"^) include both the first- and the second-order parts. 

We follow both approaches here. The PPM approach has the advantage that no nonlinear equations must be 
solved. All problems of convergence of a nonlinear solution are therefore avoided. The exact PM approach has the 
advantage that it does require a nonlinear solution and that it will help to build the computational infrastructure for 
full general relativity. It will be shown below that, in principle, the technical step is surprisingly small from exact 
PM field equations to those of full general relativity. However, there are several important conceptual difficulties that 
must be addressed in order to develop a fully general- relativistic periodic standing-wave approximation. These issues 
will be discussed at length in a forthcoming paper 0] by one of us (CB). 

The present paper focuses on the technical and computational challenges of the transition to general relativity, but 
it cannot ignore conceptual issues entirely. This is because two central conceptual difficulties of that transition, the 
problems of satisfying the gauge condition and of computing the source motion, arise already at second order in our 
post-Minkowski expansion. In fact, these two problems are closely intertwined and, in the second-order problem at 
least, can be addressed completely at the analytical level before a single line of code is written. We therefore limit 
our discussion here to just those aspects of gauge conditions and source motion that are relevant to the second-order 
theory. 

The field equations (|4|) arise by assuming the harmonic gauge condition in the Einstein equation. This eliminates 
several troublesome terms and yields an equation that, in its linearized form at least, can be solved using standard 
Green-function techniques. However, a solution of the gauge- fixed field equation will only solve the original Einstein 
equation if it happens to satisfy the harmonic gauge condition. This is not guaranteed. One can show [H that a 
necessary condition for a solution of Eq. ([4]) to satisfy the harmonic gauge condition is that the binary point sources 
generating the field should satisfy conservation of energy-momentum, VaT"'' = 0, to second order in perturbation 
theory. In general relativity, of course, this condition naturally incorporates the interaction of each source with the 
other because the derivative operator Va depends on the gravitational field. As a result, unlike Newtonian gravity. 
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conservation of energy in general relativity dictates the dynamics of the sources. That is, Va T""^ = implies a 
relativistic generalization of Kepler's law, which in Newtonian theory can be written in the form 

GM _ 

4ac2 " c2 ■ ^ ' 

This result can be considered to be the lowest-order approximation to a general-relativistically correct formula for 
some carefully defined mass as a function of a and Q. Note, however, that to derive even the Newtonian equation of 
motion, it was necessary to calculate fields to second order in M. 

To simplify the computations here we make an additional assumption, though one that is appropriate to the 
binary configurations to which the PSW approximation applies: along with our post-Minkowski approximaton on 
field strengths, we make a post-Newtonian-like expansion in orbital velocity. In particular we consider the j (? on 
the right-hand side of Eq. ([5]) to be small, and we keep only the dominant terms in this parameter. This simplification 
will apply only to the inner boundary conditions we use and to some of the details of the equations used for numerical 
computation. Aside from these points, the methods developed here are independent of this low- velocity approximation. 

The presentpaper will make frequent reference to the infrastructure built up in previous papers, in particular in 
Refs. 7] and 8]. The details will not be repeated here, but it will be useful at the outset to point out the connection 
to and differences from the Minkowski background theories of previous papers and the approach in the present paper. 
It will also be important to define several coordinate systems closely related to those of Refs. 0, [!]• We will assume 
that there exist coordinates t,x,y,z that cover the region of the manifold outside the sources, the region in which 
we shall do computation. For convenience, we shall call this system our Minkowski-like coordinate system, although 
in the present paper we do not really consider a Minkowski background. In terms of these t, x, y, z coordinates, we 
assume that the Killing vector of our helical symmetry has the form 

dt + ^ {xdy - yd^) (6) 

where 17 is a constant. 

We introduce three other coordinate systems that are also related to those in Ref. [1]. In that paper they were 
alternative coordinates for the Minkowski background. Here they are defined only as specific transformations of the 
coordinates i, x, y, z. The first of these is the system t — t,x,y,'z = z, defined by 

X = X cos fit + y sin Qt y = — x sin Qt + y cos fit . (7) 

Since C^x ~ C^y = C^z = 0, the x, y, z are labels on trajectories of the helical Killing congruence and we are justified 
in calling them rotating coordinates. Another set of rotating coordinates is the cylindrical system r, z, ip defined with 
r = \/a?^r"y2~ = -^/J^ + y^ , and with tanip = y/x. It should be noted that the Killing vector of Eq. ^ can be 
written as 



dt + fl [xdy - yd^) = 9f + fld^ ^ c%. 



(8) 



where the last expression is the derivative with Lp held constant. For scalar functions the imposition of helical symmetry 
involves the replacement 



dt - fld^ = ~fl [xdy - ydoo) = -fl {xdy - yd^) 



(9) 



It is useful, as in Refs. [3] and to define yet another set of rotating Cartesian-like coordinates X^Y^Z as a 
simple renaming 



Z — X X = y Y ~ z . 
Our adapted coordinates x, 8, $ are related to the AT, F, Z system by the transformation 



Z 
X 
Y 



a?+X^ cos 28 -f ^(a4 -f 202^2 cos 26 + 



-a2 - x2 cos 28 + ^(a'^ + 2a'^x^ cos28 + x'') 



-a2 - x2 cos 28 + + 202^2 COS28 + X'') 



cos $ 
sin $ , 



(10) 

(11) 
(12) 

(13) 
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and the inverse transformation 



Z-a] +X^ + Y^ 



Z + a] +X^ + Y'^ 



1/4 



e 



1 _i / 2zvx2 + y2 

- tan — — 

2 \ Z2 - a2 ^ X2 - y2 



$ = tan^i (Y/X 



(14) 

(15) 
(16) 

The remainder of the paper is organized as foUows. In Sec. [IT] we introduce our formal expansion of Einstein's 
equations and we discuss the truncation of this system to second order. Here it is demonstrated how in our formahsm 
the change from second-order post-Minkowski field equations to those of full general relativity involves only very minor 
modifications. We then discuss the approximation of small orbital velocity that we will use in the computations. In 
this section also, we present the derivation of second-order correct inner boundary conditions. 

Section [IIII presents the formalism underlying a numerical approach. Since we can compute only "helical scalars," 
unknowns that are functions only of rotating coordinates, we show how the techniques introduced in Ref. [S] can be 
extended to the post-Minkowski equations (and to the full Einstein equations). In this section we also give the detail 
necessary for casting the computational problem in terms of the adapted coordinates that proved to be very efficient 
in earlier work[3,@|. Section HVl deals with the numerical solution of Eq. Numerical results are given along with 
a discussion of the limits of validity of the results and the importance of nonlinear contributions. Section |V] gives a 
summary of the step taken in this paper and relates it to what remains to be done. 

Throughout the paper we adhere to the conventions of the text by Misner et aLfloj. In particular, we use units in 
which c = G — 1. 



II. THE FORM OF EINSTEIN'S EQUATIONS 



A. The full equations 



We follow here the convenient formulation of Landau and Lifschitz [111 . Il2l | for the Einstein equations. This 
formulation encodes the geometric information in the densitized metric 

g"^'^/^g"/'. (17) 

From this we introduce the Landau-Lifschitz quantities 

and 



+ 9^.. 9^' S"^,A ii^".p + i (2 ff"^ - 9"^^ 9^"") (2 gxp 9-^. - 9Xa Sp^) S^^a^ 5"^. (19) 

(Our notations differ slightly from those of p^, where ^"'^^'^ is denoted Jp"^'''' and the Landau-Lifshitz pseudotensor 
jt"'^ is denoted tt"''. We prefer our notation in the second case because we associate C"'^ in a related paper ^ with 
the true, material stress-energy. We also reserve ^"'^^'^ for other purposes.) In terms of these quantities, the Einstein 
equations G^n = iirT^i, are written as 

= 1°^ - 167r/^fet^T"'^. (20) 
Putting the definitions P^ . ([TO|) into the left side of the field equation ([^0]) . and rearranging indices slightly, we find 

V ,pi> — SS U .pi' 4J SS .pi' U .pi' U iJ U .pf T^^SJ ,pU ,v U ,v ^ ,p iJ ,pU \^^) 

We take our t, x, y, z "Minkowski" coordinates to satisfy the harmonic gauge condition that 

S"^/3 = 0. (22) 
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in these coordinates. This choice greatly simplifies the field equations, which become 

+ 5^" r ^ A 5'"^ ,p + i (2 "'^ s''" - 5"'' 3^") (2 gxp g^a - gxa gp^) e^" e"'' • (23) 

To simplify this result further we collect terms on the right hand side as follows 



+ i (2g"- - g-^) (2g^^ g^^ - g,^ g^^) ] s^*.« fl^^. , (24) 

and we rewrite this as 



per ap _ 
U 15 ,pCT — 



2 5 9p<y 



5P 61 S; 6^, ~ 2 SP g-- g^x + 51 g^^ g^^ - \ g"- g^x g^^ g""" 

+ \ g'"^ g"" i^grx g^p - gr^ gxp 



S^'^.kS^'^,.. (25) 



It should be noted that the ordinary covariant metric gap and its inverse g"^^ enter (j25p only in complementary pairs. 
We define the inverse of our basic field g"'' by 



S 5/37 = 7 



so that 



1 



gap ■ 



Using this definition, we can rewrite the field equation in the form 



(26) 



(27) 



pa aP _ 
U U ,pa — 



5P 51 <5^ <5« - 2 5P tf- S0A + 5P 51 

- \ 5"" SrA S""^ + i S''" C""^ (2erA 0^^ " 5x0 Sa^)] s'^,. 5^^.- (28) 



We next define h by 



S 



.aP _ 



v/- det77(?7"''-/i"'^). 



(29) 



Note that the determinant factor is unity in the coordinates (t,x,y,z) in which we take tensor components, but 
technically is needed to define the perturbation tensor field h"^ . In terms of this new quantity, the field equation can 
be rewritten as 



S'r 51 5^^51-2 5P 5- 5^, 5^ + 5P 51 0-^^ 
\ S"" 5.A 5"" + i e"" S"" (25rA 80^ ~ 8^0 Sap)] ^"^k + ^""/^^^p. • (30) 



It is this form of the Einstein equations that we shall use through the remaining steps of our program. Equation (j30p , 
along with the definition (|^^ . are to be considered as equations for the unknown fields h"'^ . Note that so far this 
equation is exact; there has been no reference to a split into background and perturbations. 

B. The post-Minkowski truncation and our expansion 

If we consider Jf^ a perturbation of Minkowski spacetime, then the left side of Eq. (I30p . rjP'^h'^^p^, is linear in this 
perturbation while the right hand side is of second (and higher) order. 

Keeping only the linear terms gives us the equations of linearized general relativity 







(31) 
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in which we have adopted notation like that in Eqs. ([2]),([3])- This is equivalent to the usual formulation of linearized 
general relativity since it is easily shown that, to linear order, h"^ is the familiar trace-reversed metric perturbation 



Sgpaj ■ (32) 

Thus, the h"'^ defined in Eq. ^ agrees to first order with the well-established notation for linear perturbation 
calculations that was used in Ref. [1]. Note, however, that its relation to (5g^j/ at higher order in perturbation 
theory is non-linear. 

Keeping terms to first and second order gives our post-Minkowski approximation: 



where 



'-'ripXp 



□ h"" = S%^; h^'t'^^ h^'^,, + hP'^ h'^^^p, , (33) 

- iv'^ VrXV^'pV^" + 1^'"' V^" {'^VrXV^^ - Vr4,VXp)] ■ (34) 

It should be noted that the conversion of these post-Minkowski equations to the equations of full general relativity 
requires only the replacement of the 77 s by 5 s in S"^^l^. No changes need to be made in the wave operator on the left 
or the differentiated fields on the right. This means that a computer code designed to solve Eq. (|34p can very simply 
be converted to one that solves the full theory. 

As pointed out in Sec. U there are two ways in which the equations of Eq. ([55)1 and ([M)) can be approached. Here 
we describe the PPM method, the simpler standard path of solving Eq. (|3ip first for the hs correct to first order, 
then using these first-order correct hs to construct a known right hand side of Eq. p3|) . In the notation of Eq. ([3]), 
our problem becomes 

where hP"' is a solution of Eq. ([51]) . 

To proceed, we must more carefully consider just what the nature is of our approximation scheme. In the usual 
post-Minkowski theory[l^ an expansion in field strength is used, and the particle velocities are considered to be 
of zeroth perturbative order. Here we are using a different scheme which is best thought of as as an expansion in 
the source mass M . A more careful description of our approximation is that we are considering a family of helically 
symmetric solutions of the Einstein equations describing two co-orbiting "particles," each of mass M, moving opposite 
to one another on a common circular orbit of radius a, and coupled to standing waves. The parameter e = M/a is 
the parameter on which we base a small parameter expansion. Via Kepler's law, or a relativistic extension of it, the 
velocity v of our source objects is of order y^M/a . It is convenient, therefore, to consider a factor of v to represent 
a half-order in our expansion scheme. A quantity proportional to vM/a, for example, will be considered to be order 
1.5. 

In this expansion scheme h"^ is, in principle, of first-order in M/a (or second order in v). Not all "linear" components 
are of this order, however. The linearized field equations of general relativity, in the harmonic gauge, are 

nh°'f^ = -IGnT^f , (36) 

and the stress-energy component is proportional to the source mass M and, to lowest order, is independent of v. 
All other components of the stress energy are, to lowest order, proportional to one or more factors of v. Components 
T** (where j is a spatial index) are of order 1.5 in M, and components T*^ are of order 2. Thus ft,", in linearized 

general relativity, is the only component of h"^ that is driven by a first-order source, and hence is the only component 
that is actually of first order. Here, we summarize the somewhat complicated situation regarding orders in linearized 
theory and PPM equations. 

(i) For ft**, the first order fields are found from linearized gravity, Eq. (pij) or and are of order v^M/a (that is, 
first order). Equation ([35|) is used to solve for ft**, with contributions on the right only from ft**. 

(ii) For ft*' the lowest order fields are of order 1.5. To solve for these fields, we need only use linearized theory, 
i.e. , Eq. (|3T|) or (|36p. In principle, we could adapt Eq. (pS]) . with^ft*' on the left, putting on the right products of ft** 
and^ft**, both from linearized theory. We do not solve for these corrections to ft** in the present paper since they are 
higher than second order. The computation of the lowest order fields has already been described in Ref. There 
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Eq. p6)) was solved for all components of hP'^ , and it was pointed out that the procedure is inconsistent, i.e. , that 
the solutions are of different order. 

(iii) For If^ the lowest order fields are of order 2. Thus, it is inconsistent to solve Eq. ((3T|) or ([36|) for these lowest 
order fields. The consistent procedure is to use an adaption of Eq. ([55)1 with K^^ on the left, and with products of h}'^ 
on the right. In the case of Eq. ([M]) . the spatial components of the stress energy are to be included. 

We should note that our approximation scheme has some of the spirit of a post-Newtonian rather than a post- 
Minkowskian perturbation methods. But in our scheme there are waves at every level of approximation, and there 
is no c ^ oo Newtonian limit. It is, therefore, justifiable to consider our approach to be a type of post-Minkowski 
expansion. 

C. The gauge issue 

To specify a solution, we must, of course, add source terms or boundary data to Eqs. ((3T|l or (|34p . The solutions 
thereby determined must satisfy the gauge condition in Eq. ([^ or equivalently must, in principle, satisfy /i"^ = 0. 
In practice, this condition is relaxed in a post-Minkowski approximation. If we are computing the fields correct to 
order n, then the gauge condition need be satisfied only to order n, that is: 

h'^P^p = O [h^'') . (37) 

Thus in linear theory we must only have that h"-^ ,[j be of second order, and in the post-Minkowski problem of Eq. ([M]) 
we must only have the gauge condition satisfied to second order. 

For sources moving in binary orbits, this issue is apparent in linearized theory and the way in which we deal with 
it was raised in Ref. @; we review that argument here. In our present notation, in the linearized theory, the field 
equation becomes that of Eq. ([55]) . The stress energy tensor for a pair of binary point masses (the stress energy used 
in Ref. Q) will only satisfy T"^ = if the masses are at rest. For masses moving at velocity v there are components 

of T^p that are of order vM/a and of order v'^M/a, and the solution to Eq. (|36p . therefore, cannot satisfy Jf'^ = 0. 
The way in which we are to understand this is by considering the missing terms in 

= -167rr"^ -I- O (h^) . (38) 

The second-order terms on the right can be thought of as representing the gravitational forces that drive the source 
masses. The gauge condition to lowest nontrivial order will involve the divergence of the stress energy and the 
quadratic post-Minkowski terms. Satisfying that gauge condition is what gives us the relationship of v and M/a, i.e. , 
it is what gives us Kepler's law. To the order at which we are working, keeping terms only of order (h"^)'^, we will 
infer only the standard Newtonian Kepler's law. A higher order treatment would lead to relativistic corrections of 
the relationship of v and M/a. 

D. Boundary conditions 

At distances r much larger than the orbital radius (and therefore much larger than the source masses M) the 
post-Minkowskian corrections to the metric are tiny. The boundary conditions on the fields, then, are just those of 
linear theory, as described in Ref. 8] . These conditions are the usual Sommerfeld conditions, that in the Minkowski 
coordinates all fields f{t, x, y, z) far from the sources obey dtf — drf for ingoing waves, and dtf = —drf for outgoing 
waves. (Here r is the radial coordinate y^x^ + y^+z^ .) In general relativity, of course, neither ingoing nor outgoing 
helically symmetric waves are possible, but these conditions are needed in constructing the standing wave solutions 
that are possible in general relativity. 

In order to complete the specification of the fields, we must give the conditions at inner boundaries, i.e. , at small 
coordinate distances R from the location of the sources. Here "small distances" means that i?/a <C 1 where R is some 
characteristic distance from the location of the source. Within the scope of our approximation we are then looking 
for conditions near the source that correspond both to R/a <SC 1 and to second order in M/a. 

The equations for which we seek boundary conditions are those of Eqs. p3|) and (p4)l . The right-hand side of 
Eq. (j33p is known from the solution of the linearized problem. The only ambiguity then is what homogeneous solution 
of Dft,"^ = to choose. To lowest order in M/a, the answer is clear: we choose the same moving monopole solution 
as in Ref. . The choice of the homogeneous solution to second order in M/a is more subtle. In principle we could 
have a solution that corresponds to a moving dipolc. For example, we could have 



h" =p-R/R^ , 



(39) 
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where R is the distance, in a local (approximate) Lorentz frame from a source location to a field point. The dipole 
constant p would have to represent the acceleration of the local Lorentz frame since, for a spherically symmetric source, 
there is no other spatial vector related to the physics of the problem. Dimensionally this second-order term would have 
to be proportional to M'^aP . This would not make physical sense since this "acceleration term" would be independent 
of the orbital radius a. The arguments against second-order homogeneous solutions of higher multipolarity are an 
extension of this idea. It is worth emphasizing that the choice of the inner boundary condition is equivalent to the 
choice of the "particles" whose motion produces the fields. We are free to choose both first- and second-order inner 
boundary conditions to model sources with large intrinsic multipoles. Such intrinsic multipoles would be independent 
of the orbital radius a. 

The particular solution of Eq. ([55]) can be found, of course, from the known form of the right-hand side of the 
equation. From dimensional arguments alone, one would expect that in the limit of small distance R from the source, 
the second-order solution has the character M'^/R'^, M^/Ra, and so forth. We shall impose our inner boundary 
at R <^ a, so that the second-order term M^/Ra is negligibly small compared to /R?. (For a post-Minkowski 
approximation to be valid in the computational region, we must, of course, also require that M/R be small.) We 
repeat here the conditions on the choice of the parameters of the physical configuration and the inner boundary: 

M/R < 1 i?/a < 1 . (40) 

With these choices the inner boundary conditions will be independent of a, and hence the second-order correct 
boundary conditions can be found from the Schwarzschild solution expressed in coordinates T, X, Y, 2'[l^ that satisfy 
the harmonic gauge condition [T^. [l5| 

where R? = -I- + . (Note: These particle-centered coordinates are not related to the X,Y ,Z coordinates 
of Eqs ^_([T0|) -(fT6 | . ) The next step is to expand this solution in M/R, adopting the method of Johnson and Van 
Meter[lJ|. We let z'^{t) be the coordinate path of one of the source particles, where r is the particle proper time, 
and we let = dz'^/dr be the 4-velocity of the particle. 

We next consider the metric in Eq. (|4T|) to be written as y/\detg\ g°'^ — rj"'^ — h"^ , as in Eqs. (fTT]) and (|29|) . so 
that h"^ is of first and higher order in M/R. For an event at a;^, we define 

= a;^ - z^(t) (42) 

so that is a null vector, and t is the retarded time for the event, to zeroth order in M/R. Next we define 

n = -U^r" (43) 

to be the distance from the particle to the event, to zeroth order in M/R, measured in the particle comoving frame, 
and we define 

n^' ^ — -W (44) 

to be the spatial direction to the event at x^, to zeroth order, in the particle comoving frame. In terms of these 
quantities, the fields near the particles, correct to second order in M/TZ are 

/AM 7M'^\ M^nf^n" 

+ (45) 

As in Ref. [1] we choose the sources to be locally-spherical points on circular paths One of the source points 
( "particle 1" ) moves on the path 

X = acosflt y — asinQt z — 0. (46) 

and the other ("particle 2") has opposite signs in the formulas for x{t) and y{t). The instantaneous 4-velocity is 

C/° = 7 Uy = ±v"f . (47) 

Here and below the upper sign refers to particle 1, the lower to particle 2. The particle motion clearly satisfies the 
helical symmetry, a necessary condition for the "particles" to be sources of helically symmetric fields. 
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We define v to be the coordinate speed a^l of the particles, and 7 to be the associated Lorentz factor 7 = 1/ — . 
To proceed we follow the derivation of Eqs. (22)-(28) of Ref. Q to relate the retarded quantities to our Minkowski-like 
coordinates. This is done by introducing instantaneosly corotating coordinates x, y, z: 

y = y ± vt x = x a = z = (x — a)'^ + y'^ + . (48) 

For an event t,x,y,z near (i.e. , at small TZ from) a particle it was shown in Ref. that ipart, the retarded time at 
the particle, is given by 



t'part 



t =F v^'^y — 71/ + 72^2^2 , (49^ 



From this we have 



n=^JP+ -i^v^y^ , (50) 

and for the null vector of Eq. we have 

7"^ = {i - tpait, a; =F T fipart; z . (51) 

With fpart eliminated, this is 

— {±D7^y + 77?,, x^ a, 7^2/ i ^177?,, z} (52) 

and the vector of Eq. ([ii)) is 



n 



W7^i/ X =F a 7^y z 



r_i_ ; a I / a i 



With these expressions we can now explicitly evaluate the components of the inner boundary condition in Eq. (|45p : 

-^^ /4Af 3 \ 3 2 4 



.M2 3 ^(0: to ) 



-^t.^-^.t = T;^^'7^^^^ (55) 



- _ M2 (x T a? 

" ^^^^ 

= = -^7^fe^ (59) 



S,.,5.» = -^i'Ij (62) 

A'P z^ 

III. EQUATIONS FOR COMPUTATION: HELICAL SCALARS AND ADAPTED COORDINATES 

A. The equations in helical scalar form 

In Ref. Q the metric perturbation fields were considered to live on a Minkowski background. We introduced a set of 
basis vectors n = 9t, ex = i9x, Gy = 9y, = 9z that were covariantly constant in that Minkowski background. In the 
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tnn — 


nn 






tnO = 


1 






tri,±l = 


2 1 




- {e^ ± iey)n] 


to,o = 


1 


\^X^X ~t~ GyGy 


+ eze^] 


t2,0 = 


-1 

Ve 


\Gx^x ^" ^y^y 




t2,±l = 


t\ 




^ 2 ^ ~l~ ^z^y 


t2,±2 = 


\[e 


x^x ^y^y i 





present paper, and in full general relativity, it is no longer convenient to consider fields on a Minkowski background. 
It is, however, important to use the infrastructure of Ref. [1] for constructing helical scalars. To do this we will use 
the same symbols as in Ref. Q , but with a somewhat different meaning. Here these quantities are to be interpreted 
as indexed objects whose components are constant. Thus, for example n = c)t, is to have indices — 1, nt = —1 and 

= rij = 0. The index manipulation of these symbols will done as if they were components in a Minkowski basis. 

In Ref. 3 , we constructed 10, rank 2, symmetric basis tensors that were covariantly constant in the Minkowski 
background. Here it is convenient to use the same notation 

(64) 
(65) 

(66) 
(67) 
(68) 
(69) 
(70) 

Here we intepret these simply to be constant indexed objects with the numerical values they would have in a Minkowski 
basis. Thus, for example, = 1. Lastly, we define 

tnn = t„„ t„o = tjiO too = tpo t2Q = t20 (71) 

t„,±i = eT^"*t„,±i t2,±i = eT''"*t2,±i (72) 
t2,±2 = eT2*^"t2,±2. (73) 

We can express h"^ as a sum using either type of tensor-like basis: 

/i"/3 = *^t"^ = $^t"^ = ¥°"'^t^^ + ... + ¥^'-^^tf_^ , (74) 

where the label A takes any of the 10 values {nn) • • • (2, —2). 

The components t'^^ are not constant, but they have a very useful property. They behave under £^ like the 
components of a Lie-dragged tensor. This means that the scalar-like quantities constructed in Eq. ([7^ will be 
helical scalars, i.e. , they will be constant along the helical trajectories. As in Ref. these basis tensors have the 
following important property for differentiation with respect to the Minkowski coordinates 

df^tf - d,, (e-^M^)nt ^a/s^ ^ ^^^^^^ ^-.^(A)ot ^a/J „^ ^ ,^(A)f] tf , (75) 

where ij,{A) has the value for A = {nn), {nO), (00), (20), has the value ±1 for A = {n, ±1), (2, ±1) and has the value 
±2 for A= (2, ±2). 

The representation for h"^ can be substituted into Eq. and the result can be contracted with the (orthog- 
onal) basis symbols t"^. The result is a set of equations of the form 

□5"^ - 2ifi{A)n^d^^'^ + n{A)^n^^'^ = . (76) 

For A = (n, ±1), (2, ±1), (2, ±2), the fields 4'"^ are complex, and in practice we work with the real and imaginary 
parts 

= C/^ + iV^ , (77) 

of these fields. Separated into its real and imaginary parts, Eq. (ffS)) becomes 

aU"^ + 2n{A)n^d^V^ + n{A)^n^U^ = Real part of (Q^) (78) 
□y^ - 2fi{A)n^d^U'^ + ^(A)2f7V^' Imaginary part of (Q^) . (79) 

The □ operator here, as in Eq. (|76|) . is rj'^^dadp, and the helically symmetric time derivatives are implemented through 
the replacement dt —n{xdy — ydx) — —^d^, so that 

□ = a| + a? + dl - v?dl . (80) 
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In Ref. [31 the right-hand sides Q were zero in the region outside stress-energy sources. In the post-Minkowski 
approximation of the present paper, the right hands side follow from the contraction with t"^ of the right-hand side 
of the equations in ([55)1 , , expressed in terms of the representation of h°'^ given in Eq. (fft)) . The effective source 
term Q-^, then will consist of terms quadratic in derivatives of the 4^"^, and [from the last term on the right in Eq. (|33p ] 
in products of ^I''^ and its second derivatives. The form of the source terms then is 

= S'^'ic da^"" ^h*"" + ^"'sc ^a^^*^ (81) 

where the coefficients S°'''^(j, T'^^gfj, relatively simple functions of the rotating coordinates, are derived in the 
appendix. Equations (j76p and (|8ip are the field equations we solve for our post-Minkowski approximation. 

The inner boundary conditions, are those of Eqs. ([51)1 " converted to projections on the helically symmetric 
basis vectors. A straightforward computation gives 

*"°=±^^ («) 

. T (1^ + fg) W T ^ ^ (84) 



~2n /4M 7M2\ w2 2 m2 (x T a)^ + 7"*^^ - 222 



* ^ 



Here and below we omit equations for A = {n, —1), (2, —1), and (2, —2), since the quantities carrying these labels are 
(up to a sign) the complex conjugates of the the quantities with A = (nl), (21), and (22). 

In Sec. H] it was explained that our expansion is really based on a family of solutions of Einstein's equations with 
varying e = M/a, with M the mass of a source, and a the radius of the orbits. We noted that the orbital velocity, 
of order y/M/a , can be used to keep track of orders of terms, and we noted that of the components h°'^ only /i" is 
truly first order. 

That argument was based on considerations of stress-energy sources in linearized gravity. In the present paper, 
we represent the effect of the source objects through inner boundary conditions, rather than through explicit stress 
energy sources. In Eqs. we can see explicitly the orders of those inner boundary conditions. Again, of the 

components of h"^, only h** has a first-order piece to its inner boundary conditions, a piece that goes as {M/a)v^ . 
(Note: The factor M/TZ should be viewed as the order parameter M/a, divided by the dimensionless distance TZ/a 
that is the limit appropriate to "near-source" inner boundary conditions.) The inner boundary condition for ft,** also 
has a second-order piece M"^ /a? and a piece {M"^ /a?)v'^ that is third order in M/a, or sixth order in v. By similar 
counting, we see that the lowest order term for h}^ \s ^ , while /i*^ and h}^ are ~ w^, with all other components 
w"*. For the description of the field in terms of helical scalars, Eqs. ([55)1 - (|10ip tell us that is of first (and 

higher) order in M/a or of order , while is of order and all other helical scalars are at least of order 

We summarize here the leading orders of all 
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A consequence of this is that for our second-order post-Minkowski solutions only the ^I/^"") terms need to be kept on 
the right in Eq. ([5T|) . so that those source- like terms can be simplified to the order w"* expression 

The details of these terms are given in the appendix. What is of greatest importance is the order of these driving 
terms: 

Q(-)^„4 Q("0),g(»l) g(00)^g(20)^g(21)^g(22) _^4_ (9^) 

It follows that the computation of the second-order vpC"") requires a solution of a nonlinear problem. To lowest order 
(w^) computation of is linear, and does not require a driving term. For the other fields ( ^-^^o), 

\[/(20)^V[r(22)'j^ the driving term is of the same order as the lowest order field. In this case, the solution for the lowest 
order field requires inclusion of the driving term. The resulting problem is not nonlinear, however, since the driving 
term involves not the field being computed, but rather the first-order field 



B. Adapted coordinates 



The corotating coordinates x, y, z are not well suited to describing the fields and sources of the rotating binary. As 
in Refs. 0, Q we introduce a set of "adapted coordinates," Xi ^ transformations of x, y, z that are better suited 
to encoding information about the sources and fields, especially in combination with the truncation of a multipole 
analysis, as laid out in Refs. 0, 

The introduction of adapted coordinates imposes two changes in the details of Eqs. ([75)) . ([5T|) . First the operators 



□ and d^p take the form 



- ^ee 902 + 



g2^ 

+ ^xjr^ 



d_ 



^= r 



96 



dx 



(93) 



where the coefficients B, F are known functions of x, O, $ given explicitly in the appendix of Ref. @]. The second 
change needed is that the derivatives in the expressions for the must be converted to derivatives with respect to 
the adapted coordinates. Since the rotating coordinates are relatively simple functions of the adapted coordinates, 
this change is straightforward; the details are given in the appendix. 

Though each step of the transformation to adapted coordinates involves elementary functions, the full set of steps 
that must be taken in transforming, and in imposing the eigenspectral method and multipole filteringf?'! becomes 
exceedingly tedious and prone to error. For this reason, the transformations, eigenspectral method, and the generation 
of the final finite difference equations for computation have been implemented as symbolic manipulation with Maple. 

With the helical symmetry of Eq. ©, and with the fact that far from the sources x ~* ''i the Sommerfeld outer 
boundary conditions, dt"^^ — ±3^"^^, becomes 



YX 



dx 



(94) 



where -I- and - correspond respectively to outgoing and ingoing conditions. 
The inner boundary conditions of Eqs. (|82p - (|88|) become 



/4M 

[iz 



7M' 



M2 z;2^4 



7^2 7^2 4a2 



sin^ (29) cos^ $ 



(95) 



(96) 
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= 





[AM 


7M2\ 









TV 7^2 4a2 



sin(29) cos $ [sgn[cos 9] cos 28 - 27^ sin 26 cos $] 



= + — ^ -4^ 1= — At (1 + 7^ - 1) sin^ 26 cos^ 2$) 



^20 ^ 



/4A// 



7M2\ ^272 



V 7^ 7^2 ; ^/6"7^4 4a2 

M2 



21 



22 



27^4 4a2 

Here, 7?., in terms of adapted coordinates, is given by 



sgn[cos 6] (97) 
(98) 
(99) 
(100) 

cos2 26 - 7* sin2 26 cos^ $ - sgn[cos 6]2i72 cos 26 sin 26 cos $) . (101) 



M2 



^ (cos^ 26 + 7^ sin^ 26 cos^ $ - 2 sin^ 26 sin-" $) 



^. ^ sin26sin<i> (cos26sgn[cos6] - i7^sin26cos$) 
7c^ 4a^ 

/4M 7M2\ ^2^2 ^2 ^4 



7^2 



4a2 



[1 +72t;2 sin2 26cos2 $] 



(102) 



In principle, this completes the specification of the problem to be computed in adapted coordinates. We summarize 
that problem here, (i) Our field equations are those of Eq. for A = (nn), (nO), (00), (20) and Eqs. ^ 

for A = (n,±l), (2, ±1), (2, ±2). In these equations, □ and are given in adapted coordinates by Eqs. (|92|l . ((93l) . 
(ii) The source terms Q"* are given, schematically, by Eq. ([M)) . with indices a,b taken as adapted coordinate labels. 
{Hi) The outer boundary condition is given by Eq. ([M]) . (iv) The inner boundary conditions on the fields are 
given by Eqs. ^ - pOTj) . 



IV. NUMERICAL METHODS AND RESULTS 

It was explained at the end of Sec. IIII Al that we work in the approximation of small orbital velocity v. Of the fields 
"I*^, only requires the solution of a nonlinear problem in that small v limit. For this reason we shall emphasize, 

in the presentation of the results, those for vl/(""). 

The nonlinear problem is that defined, for A = (nn), in Eqs. (|76|) and (|90p . More specifically, it is 

□$("") ^ (^:^rf^ + in'^n'') + n^n" ¥Z''^ . (103) 

The orbital velocity v is related to the mass M of cither of the orbiting masses by 

^2 = M/4a , (104) 

and to linear (i.e. , lowest) order vI/("") ~ v"^. Here we compute correct to second order, i.e. , to order w*. 

Inner boundary conditions are imposed on an approximately spherical small surface at x = Xmin- The fact that we 
are using a second-order post-Minkowski approximation puts a significant restriction on the value of Xmin- From the 
definition of x in Eq. (fH)) (see also [7, 8]) it follows that at a small distance R from one of the source objects 



X w V2aR . (105) 

From this and Eq. (jl04p we have that the maximum field strength for the domain of computation, the field strength 
at the inner boundary, is 



M 



-^v^^^^n^. (106) 

Amin Amm 



For the second-order post-Minkowksi approximation to be justified, this measure of field strength must be significantly 
less than unity, so for a given choice of we must in principle choose the location of the inner boundary to satisfy 



Xmin/a > 2\/2 an = Xmin/a > 2V2 v . (107) 
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In the results to be presented we will vary and Xmin to achieve different values of M/R\^^^. This will allow us to 
compare the computed errors due to the post-Minkowski truncation with the expected errors. 

An estimate of errors will be possible because we will present the results of three different approaches to the 
computation of vl/(""): 

(i) A computation using linearized general relativity, precisely the computation done in Ref . |S] . This computation 
is done by solving □']/("") — for the lowest order ( — AM /TV) boundary conditions in Eq. ((95|) . imposed at 
Xmin- This should give correct to first order in M/R\^^^. 

(ii) A computation using the perturbative post-Minkowski approach, as indicated in Eq. ([3]). This approach starts 
with the hnearized computation of (i) above to find ^F^""). That first order result is then used in Eq. ()103|) in the 
form 



(nn) 



{^^r^^'^ + \n''rf^ + n''n'^*("")f . (108) 



This equation is solved for »!'(""), using only the second-order part of the boundary condition in Eq. I|95p . (That is, 
the linearized part of the boundary condition, i.e. , the part of the boundary condition used for ^("■"), is subtracted.) 
The numerical result involves the same "eigenspectral" method as was used in Ref. [1] , with the addition of the driving 
terms. The final result is 

This result is expected to be correct to second order in M/R\^^^. 

(iii) A computation using the exact post-Minkowski approach. Here Eq. (jl03|) is solved as a nonlinear equation, 
and the full boundary condition in Eq. (I95|) is used. The equation is solved by a multidimensional Newton-Raphson 
scheme similar to that used for the nonlinear models in Ref. Q . The result of this computation is expected to be 
-^{nn) correct to second order in M/R\^^^. Since this exact PM approach, and the PPM approach described above, 
are both second order computations, we expect the difference between the two methods to be third order in M / R\^^^y^. 

The results for several models are shown in Figs. [1] - [31 Since the monopole moment is much larger than the 
radiation fields, these figures show the computed results with the monople moment subtracted. The sharp feature 
near Xmin shows that the inherent quadrupole moment imposed by the inner boundary conditions is immediately 
overwhelmed by the quadrupole moment due to the binary configuration. This near irrelevancy of the inherent source 
structure has been discussed in detail in Sec. V, and Appendix B, of Ref. [3]. 

Each of Figs.[T]~-[3]corresponds to a choice of binary velocity v — aVl, and of the inner boundciry paicinieter Xmin- 

All 

computations were done with a 16 x 32 x 16000 grid respectively in 8, $, x, and with six (discrete) spherical harmonics. 
The Newton-Raphson iteration of the exact PPM computations were iterated 20 times, although convergence was 
typically achieved after only three or four iterations. Each plot shows the results of three different computations: those 
for linearized theory, perturbative post-Minkowski, and exact post-Minkowski. Figures [1] and [21 for orbital velocity 
v/c = 0.075, shows the effect of varying the choice of Xmin/a through the set of values 0.3,0.45,0.6,0.7,0.8. The 
good agreement of all three (exact, perturbative, and linearized) computations for Xmin = 0.6, 0.7, and the excellent 
agreement for Xmin = 0.8 is an indication that for these choices the inner boundary is large enough that the field 
strength is small at the inner boundary, and the post-Minkowski approach is justified. For Xmin = 0.3 there are 
significant disagreements among the results of the three approaches. For this case, the crude estimate in Eq. (|106[) 
gives a value of 0.5 for MjR^ too large for a weak field approximation to be reliable. 

Figure [H shows a detail in the wave zone, for the five v/c — 0.075 exact computations. The results shows that 
the strong field error for Xmin/a = 0.3 significantly reduces the wave strength. For the larger values of Xmin/a 
the differences in wave strength are small, and have no pattern, indicating that the error is being dominated by 
sources other than the strong field at the inner boundary. Figure [3| shows computational results for the larger orbital 
velocities 0.1 and 0.15. In these examples we see further evidence that higher orbital velocity, and smaller Xmin 
leads to disagreement of the three (exact, perturbative, linear) computations, suggesting a violation of the weak-field 
requirement. 

Table [J gives a crude analysis of the correlation of the field strength at the inner boundary and the agreement of 
the three computations for a given model. The parameters of the model are given in columns two and three with the 
fourth column giving the estimate of the field strength at the inner boundary, according to the criterion in Eq. (|106p . 
The remaining columns give various indicators of agreement of the linearized, PPM, and Exact PM computations. 

In the near zone the comparison is done by using differences in the value of the first maximum of \I>("") near the 
inner boundary; this is the maximum of each of the curves of Figs. [T] and [3] (The value at the inner boundary itself is 
fixed by the boundary conditions. At the inner boundary, therefore, the results for PPM and Exact PM must be the 
same; the difference between these two and the linearized computation would only reveal the second order difference 



15 




1x10 ' 



-1x10" 



-2x10 ■ 



-3x10 ■ 




60 



FIG. 2: Comparison of exact post-Minkowski computations of vI/(""'^ for aQ = 0.075, for five different locations of the inner 
boundary Xmin/a. Each curve is labeled with the values of Xmin/a (0.3, 0.45, 0.6, 0.7 0.8) for which it was computed. 



in the bounary values used.) The columns labeled "Lin v. PPM" gives the fractional difference of the computed value 
of this maximum for the linearized and the PPM computation. The column "PPM v. Exact" does the same for the 
two second-order post-Minkowski computations. The following two columns test the hypothesis that the difference 
between the linearized and the PPM results are second-order in the boundary-value field strength, and that the PPM 
vs. Exact PM results are third order. If those order estimates were accurate, the numerical values in columns seven 
and eight would be expected to have little variation. Columns nine through twelve give the same indicators as those 
in five through eight, but now for the wave zone. In the wave zone the criterion for agreement is taken to be the 
maximum of the first wave, e.g. , at x ~ 28 in the plot for the afl = 0.10, Xmin = 0-80 model in Fig. [31 

The values in Table [H for the near zone, show weak evidence for the expected effects of the field strength at the 
inner boundary The values in the wave zone show less evidence and, along with Fig. [2l suggest that the error in the 
wave zone tends not to be dominated by the location of Xmin- 

The lack of clear evidence of field-strength effects is not a complete surprise. For large values of M/R, the field 
strength for these models, we should expect the post-Minkowski approximation to be too crude. The disagreements 
on the order of 100% are in accord with this, and an analysis based on orders of a small parameter should fail. For 
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FIG. 3: Comparison of results for linearized, perturbative post-Minkowski, and exact post-Minkowski computations for higher 
velocities. 
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Lin/PPM 


Lin/Exact 
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model 


n 


Xmin 


M/R 


PPM 


Exact 


x{R/Mf 


X {R/Mf 


PPM 


Exact 


x(i?/M)2 


x{R/Mf 


la 


0.075 


0.3 


0.5 


0.61 


0.67 


2.4 


5.3 


.031 


0.85 


0.12 


6.8 


lb 


0.075 


0.45 


0.222.. 


0.34 


0.10 


6.9 


9.1 


.03 


0.18 


0.61 


16.4 


Ic 


0.075 


0.6 


0.125 


0.24 


0.006 


15 


2.9 


.032 


0.64 


2.0 


33 


Id 


0.075 


0.7 


0.0918 


.20 


0.0016 


24 


2.1 


.03 


0.03 


3.6 


39 


le 


0.075 


0.8 


0.0703 


0.17 


0.006 


34 


18 


.029 


0.020 


6 


58 


II 


0.100 


0.8 


0.125 


0.31 


0.019 


20 


9.7 


.067 


0.063 


4 


32 


III 


0.150 


0.6 


0.5 


1.0 


0.36 


4 


2.9 


.23 


1.0 


9 


8 



TABLE I: Agreement of the computations and field strength at the inner boundary. See text for details. 



significantly smaller values of M/R the error induced by the location of Xn 
dominate. 



are small, and other sources of error 



SUMMARY AND CONCLUSIONS 



In this paper we present methods and results for the post-Minkowski level of computations within the Periodic 
Standing Wave approximation. More specifically, we present the theoretical background for the PM computations 
along with two distinct ways of computationally treating the second-order PM equations: (i) the perturbative PM 
method in which first-order fields are initially computed and used as sources in the linear equations for the second 
order fields, and (ii) the exact PM method in which the second-order PM equations are solved as if they constituted 
an exact theory. 

Although the present paper deals with the second-order PM approximation, the full equations of general relativity, 
in the harmonic gauge, were given in a form easily adaptable to helical symmetry. In fact, an important point made 
in the paper is that the computational structure of the exact PM problem differs very little from that of full general 
relativity. 

In the development of the infrastructure for for the problem, we showed the utitility of the formalism of the helically 
symmetric complex field projections introduced earlier for linearized gravity theoryQ. That formalism serves 
well not only for the second-order PM work, but also for full general relativity. 

The details of the mathematical description of fields and motions clarify the difference between the PM approx- 
imation for a binary system and a post-Newtonian approximation for that system. At a characterisitic distance R 
from one of the sources, the gravitational field strength is of order M/R. Our PM approximation demands that 
M/R <C 1, since it was the very smallness of this term that justifies omitting higher order terms in the PM approxi- 
mation. The PM approximation then, in a sense, is an approximation focused on the near-source fields. By contrast. 
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the post-Newtonian approximation is one in which the source velocity is small compared to c, or M/a ^1. In a 
very relativistic black hole binary, one in which the two holes are almost in contact, there is no significant difference 
between M/R and M/a. But the PSW approximation makes no sense unless a/R is an order or magnitude or more. 

Thus, in our equations, we can distinguish the approximations associated with field strength and with small velocity, 
and we find it useful to do this. The limitation of field strength, that associated with the PM approximation, is 
determined by the value of M/R at the inner boundary. The parameter for the post-Newtonian approximation, by 
contrast, is simply (ai7/c)^. In our computations we consider values of M/R that are not exceedingly small, such as 
the M/R = 0.5 models of Table H] We may, at the same time, keep terms that are second order in v'^ /c^ ^ M/a, 
though these terms may be orders of magnitude smaller than 0.5. The choice we have made in this paper is to keep 
terms to second order both in M/R and in M/a. These distinctions will become even more important in the PSW 
approximation for full general relativity. 

An important aspect of the present paper is that it shows that there are no insurmountable computational difh- 
culties in computing the post-Minkowski PSW fields. This more-or-less guarantees that there will be no significant 
computational difficulties in the PSW problem in full general relativity, as long as we choose the inner boundary 
sufRciently large. The computational challenge for the full general relativity problem will be in dealing with the 
strong fields in the case that Xmin is chosen small enough to give a good representation of conditions very near a black 
hole. The full general relativity problem will, of course, also entail interesting issues of interpretation. 

Lastly, the fact that wc can now compute second-order PM fields means that, in principle, our results can be used 
as trial initial conditions for numerical evolution codes. In practice, our results are limited to the region x > Xmim 
and evolution codes require data also for x > Xmin- The obvious first attempt at a remedy to this is to glue a pure 
Schwarzschild puncture into the x > Xmin region. 
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APPENDIX A: DETAILS OF THE FIRST ORDER DRIVING TERM 



Here we give the details of the terms on the right in Eq. ([50)1 . We start with the Einstein equations truncated to 
second order, as given in Eq. (I34|) . and we write this equation as 

□ h'^^ = Sf^^l Ji^*,, /i^^, + hP'^ h-f'^p„ . (Al) 

As argued in Sec. Ill Bl only the /i'* terms have a nonzero first-order part, so only these need appear on the right, and 
the driving term on the right in Eq. (jAl[) simplifies to 

□ /i"'3 = S°'f^'"' h'\^ + (5f /i" h'\u , (A2) 

where 

gafiKy ^af3KU (A3) 

We next substitute h"'^ = ^-^t"^ , from Eq. on the left in Eq. (gH and use the fact that /i" = to get 

□ (^-^t"^) = t^^D^-^ = s''^"'' + ^!r"' • (A4) 

We can now use the fact that the basis tensors t"^ have the following orthogonality property under contraction 
and complex conjugation 

(t"^)*ta/3i3=<5ABe(^), (A5) 



where 



g(^) ^ J -1 if ^ = M), (nl), (n, -1) ^^g^ 
I -1-1 otherwise . 
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Using this, we contract Eq. (|A4p with {tap a)* to get 



(A7) 



Multiplying by e*^^'^)^* puts this in the form of Eq. ([7S|) . and we conclude 

This result gives in terms of derivatives of the fields with respect to the x, y, z Minkowski-like fundamental 
coordinates. What is needed is the form indicated in Eq. ([5T|l : derivatives of the fields with respect to the 
adapted coordinates. The change to this form requires two transformations. First, the derivatives with respect to 
the four Minkowski-like coordinates t,x,y,z must be changed to derivatives with respect to the rotating coordinates 
t = t,x,y,z. In doing this, helical symmetry is imposed by using Eq. ^ to replace dt- The fact that Q'^ has been 
constructed to be a helical scalar guarantees that there will be no explicit time dependence in Eq. (jA8[) : Q-^ will be a 
function only of rotating coordinates. How this comes about is related to the fact that {tap a)* = e~*''*-"^-'^*(tQ^ a)* 
so that 



,tt 



= e{A) (ta0 a) [S'^^'^'' Vl,(^-) + SA^I^nn) 

As an explicit example, here we evaluate Q("i). The general expression in Eq. (jA9p becomes 

Q^"'^ = („!))* 5"^'^'' $|r' ^!r^ ■ 



(A9) 



(AlO) 



[Recall that \I>("") = \]/("").] It is simplest to use the t,x,y,z basis for evaluation. From the definition in Eqs. ([5^ 
we have that t*-(^^-^^^ = l/\/2" and t(j|(„i) = -i/V^ so that 



g("i) = _ 



4V2 



(nn) 



5n 



4V2 



(All) 



where the time derivative of <!'("") has been replaced following the helical prescription of Eq. ^ 

It is convenient next to reexpress the results in terms of the rotating coordinates X, Y , Z in which the Z axis is 
aligned with the sources, rather than the rotational axis. For our example, this gives 



g(»i) ^ 



5n 



4V2 



■ X ,z 



.z ,x 



(A12) 



Finally we can convert the derivatives with respect to X^Y ^Z to derivatives with respect to the adapted coordinates 
X, 6, 3* by using the following relationships (see Refs. [zlilll)- With the definitions 



= ^202 + 2x2 cos2e + 2Q Sm = \/-2a^-2x^ cos2e + 2Q 

the partial derivatives for transforming to adapted coordinates take the fairly simple form 

dx cos $ (a^ + Q) 



dX 



2x3 



(A13) 



(A14) 



dx _ ■5'„, sin $(a2 + Q) 
dY ~ 2^ 



(A15) 



dx ^Sp^ 
dZ 



Q) 



2x3 



(A16) 



99 _ S'pCos$(-a2 + g) 

dx ~ v 

ae _ S'pSin4>(-a2 + Q) 
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(A17) 

(A18) 
(A19) 
(A20) 
(A21) 



5$ 



ay 




99 


5„ (a2 + Q) 




2x4 




2sin$ 


ax ' 




a$ 


2cos$ 


9r 


'S'm 



The first factor in our exampie in Eq. (jA12p has already appeared in Eq. 



,x ,z 



The r coefficients, given in the appendix of Ref. [3], are 



V ax dz) X 
~de _ ~de\ ^ cos<i>(a^cos(2e) + x") ,^25) 

dX dz) 

. ( - X^U - f^in -i>sin(29) 

* ax dZj -a2 -x2cos29 + Q 
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